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ON THE THEORY OF EVOLUTION •> . 
OF COMPLETELY MIXED STARS 
' WITH MASS LOSS 



Richard Stothera. 
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The consequences of substituting the !ft-L*%i (mass- 
luminosity-mean molecular weight) relation of stellar 
structure theory in place of the "empirical" Ift-L and 
relations are investigated for the Russian theory of 
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evolution of completely mixed stars with mass loss. ^The 
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’ mass loss rate becomes a free parameter* but stars still 
evolve along the main sequence. It is found that by 
theoretically evolving the initial luminosity function, 
the observed luminosity functions of early-type galactic 


clusters may be reasonably well reproduced with a variety 
of mass loss rates and ages. Hence the luminosity function 
per se does not give definite evidence for the mode of 
stellar evolution (homogeneous or inhomogeneous) * rate 
of mass loss# or "age of a cluster. 


1» Introduction .—- According to a theory developed 


primarily by the Russian astronomers, the upper main 
sequence on the H-R diagram is a locus of stars under- 
going evolution with complete mixing and substantial mass 
loss. Inhomogeneous evolution would take place in a few, 

slowly rotating massive stars (to explain the observed 

\ 

s 

Trumpler turnoff in open clusters) and in stars evolved 
to a sufficiently low mass so that rotation and mass loss 
are ineffective (to explain the subgiant and giant branches 
in globular clusters). However, the majority of stars 
observed today will have evolved down the main sequence 
from some initial mass lying within the range of an as- 

/ 

sumed universal initial mass function (see, e.g., Masevich 
& Parenago 1950, Idlis 1957, Masevich 1959, Fesenkov & Idlis 
1959) . 

The rate of mass, loss has been taken to be Fesenkov 's 
(1949) law 

- - k L. (1) 

dt 

and the mass-luminosity relation to be that derived from 


observations , 


From the rate of hydrogen depletion in a completely mixed 
star and a relation between the mass and mean molecular 
weight, Idlis (1957) was able to derive the value of the 
constant k. Using this value of k and integrating equa- 
tion (1)# Usher (1963) evolved Idlis' s initial luminosity 
function and found disagreement with the observed luminosity 
functions of two typical galactic clusters. In particular, 
an unobserved accumulation of stars was predicted by the 
theory. Roberts (Henyey 1960) had earlier come to a 
similar conclusion. 

However, the semi -empirical basis on which the rela- 
tion between mass and mean molecular weight depends is 
highly uncertain (cf . Section 2) . If we drop this ex- 
plicit relation and, instead, allow k to be a free 

I ' * 

parameter, stars will still evolve along the main sequence 

(cf. Section 3). In this paper we shall investigate the 

consequences of evolution in the unconstrained scheme. 

2. m-|i relation .rr.Idlis -(1957) ~USed a semi-empirical 

relation between 1ft and ii to obtain the mass-loss coefficient 

225 

k. Although this relation, ~ * (Severny 1954), was 


obtained with the help of published abundance analyses 

■ \ 

of three stars (including the sun) , the evidence is 
generally held to be inconclusive. For instance, strttmgren 

(1963) finds adequate fits to the observed main sequence 
using the same "normal " chemical composition in models 

A 

of B, A. and F stars; in particular, the position of 
Sifius A is predicted with X = 0.70. From available ob- 
servational data, Ezer and Cameron (1964) find a value 
X ** 0.74 for the sun; their use of this value gives 
reasonable results for a model of the present sun constructed 
at' the end of a series of. pre-main sequence contraction 
models. However, observations of B stars yield X = 0.60 
(Aller 1961), i.e. an increase of mean molecular weight 
with increasing mass. This is presumably due to the age 
difference between the low and high mass stars, and reflects 
the chemical evolution of the Galaxy. 

,»■ • Ideally, the to-p relation should be determined for 
w • 

a single cluster, wherein the stars are roughly coeval. 

If |i is constant and the cluster luminosity function is 
not the Idlis function, then evolution does not proceed 
with complete mixing. 

3: Evolution along the main sequence .— A spread in 

luminosity for stars of the same spectral type is observed 

* I , 

on the upper main sequence. Some of the spread is due to 

observational errors, such as uncertain extinction cor- 
rections in regions of patchy absorption. From spectroscopic 
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evidence# a high percentage of binaries occurs among 
early- type stars# causing spuriously high luminosities. 
Different rotations may also cause luminosity differences. 
Finally# evolutionary effects will result in departures 
from the initial main sequence. 

Since# in addition, our knowledge of stellar masses 
is very incomplete (Schwarzschild 1958) # the specification 
of unique to-L and relations seems to be unwarranted. 
Now exact calculations show that# for a wide range of III 
and X, homogeneous main sequence models fall close to 
the observed main sequence (BlaCkler 1958# Iben 1963). 
Furthermore# Gamow (193.8)- has -shown- that, without mass 
loss# completely mixed stars will evolve up the main 
sequence until the hydrogen has been practically ex- 
hausted* Therefore we shall adopt a fti-L-X relation as 
being more realistic for evolutionary calculations. 

The form of this relation must be gotten from theory. 
Eddington's (1926) mass -luminosity law is strictly ap- 
plicable only to stars in which nL(r)/hi(r) is a constant; 
in theory, therefore# only to gravitationally contracting 
stars of extremely high mass (and hence high radiation 
pressure) • In this case it will be given by 
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L 


m 

1 + x # 


■ (3) 


where electron scattering dominates the opacity. This 
may be considered as a limiting case. For ordinary massive 
stars with the same opacity source, the dimensionless en- 
velope parameter C (Schwarzschild & H&rm 1958) yields 


L « 


4 

-U 

1 + X 



(4) 


For stars in which bound-free absorption processes also 
occur, combination of Kushwaha's (1957) parameters A and 
C yields 



8 

_U 

1 + X 



(5) 


4 

The empirical result that L ~ to along most of the main 
sequence is confirming evidence of the validity of these 
essentially dimensional arguments. 

Comparison of the results for detailed models of 
homogeneous stars (derived from a variety of sources) 
may also be used to determine the mass and chemical compo 
sition dependences, whenever the range in X for a given 
mass' is large enough. The adopted dependences are given 
N in Table 1, where X and v are defined by 



I 


l ~ (l + x)“ v m\ 


-1 4 

It may be noted that p ~ (1 + X) * for small values of 
Z. The bolometric magnitude listed for each mass in 
Table 1 refers to an (initial) hydrogen abundance of 

about 0.7 (Henyey, LeLevier, & Levee 1959). 

\ 

4. Basic Equations . —The rate of hydrogen depletion 
in a completely mixed star is given by 


atx.mx.. It _ 


where E - 6.0 x 10 erg/gm and § is the rate of mass loss 

1 

in the form of corpuscular radiation. The full rate of 
mass loss is 


dt 5 c^ ' 


If § becomes as small as L/c (the loss due to radiant 

energy), then d Ih/dt may, as usual, be neglected on an 

2 

evolutionary time scale. Therefore we omit L/c entirely 
and obtain from equations (7) and (8) 


i3X „ _ It- 
dt EIU* 


7 


The form of § will be assumed to have a power-dependence 


on Li 


dt 


- k L 


i 


( 10 ) 


Finally, 

* 

as 


the mass -luminosity law may be written generally 


I 

! . 



l c u + x o ; v Bi 0 J : 


(id 


where a zero subscript refers to the initial epoch. 

Equations (9) , (10) , and (11) , together with the 

initial conditions 


at t - 0: m = to Q , L =. L q , X = X q , (12) 


determine the evolution of the star completely. 

5. Solutions .—Let us introduce the following non- 
dimensional variables: 


m 


m ** r , l = 


_ 1“ a. _ _ 

L ' X 1 + X ' T 
o o 


(i + x )e m 

o o 


(13) 


Then the basic equations become 


^m a . „ ,i . _ -v \ n ,s 

j _ ° > j, a K v i v a x m t (14) 


dT 


m dT 
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with K * k E(1 + X )L l “ 1 and 

o o 


at t * Os m » 1, l m i, x - 1. 


( 13 ) 


;■ I ■■ . / ,, , 

f ' ■■ : • 


In general^ equations (14) may be integrated 
analytically for m, although a simple quadrature is in 

some cases necessary to obtain t. We shall here write 

\ : 

\ 

down the general solutions. 


Case t ji 1 


l “ x“ v m* 


m • O + K {* v( 1 ' ,)+ 1 -i}] 1/ui ’‘ > i 


(16) 

(17) 




Cass. i * 1 


* - e l K(x ■ 11 


(19) 


m *» e 


K(x - 1) 


( 20 ) 


T = 


V + 1 


(1 - x v + b 


(X * 1) (21) 


or 




! • 


L (v-i 


i»0 


i)l C(X-l)K) 


1 -’ r x v ' i eO.-DKd-*) 


m[: 


(X>1V ( 22 ) 
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for integral values of v* 

If no mass loss takes place (k =» 0) , 


l * x v , m = 1, t 


tr C 1 - * v+1 > < 23 > 


This time scale is the same as for L ~ Rl (X = 1) , because 

% 

\ 

in both cases v the mean rate of energy generation, L/to, 
depends only on X. 

The case considered by Usher (1963) corresponds for- 

* r • ■ ■■ .. -V. ■ ■■ • 

mally to the neglect of changing chemical composition 

(v » 0) , 


x - 1, t - t - 1 


y ( 24 ) 


since, under the assumption of equation (2) , the mass 
dependence of X does not affect t and need not be specified 
except to determine k. 

If the mass loss is extensive enough, the luminosity 
will initially decrease. The minimum, value of K for which 
the luminosity never increases is 


U -1 

X x ’ 


Since x a 0.5, we have that K > 2v/X for L never to rise 
during hydrogen-burning. However, if K < v/Xx, L will 


10 - 


return to its initial value when x is given by 


tn x ■ * K(x - 1) . 
v 


(26) 


6. Initial luminosity function.-— We have 'adopted 
Idlis's (1957) initial luminosity function Y^^# which 

exhibits a large maximum near *» 0 and a smaller 

\ 

maximum near'Mj^ ■ - 4.5. This function agrees 

reasonably well with functions derived by Masevich (1956) 

and by Salpeter (1955). in the treatment which Idlis 

gives them. However. Usher (1963) found that Sandage's 

(1957a) recomputation of the Salpeter function predicted 

too few of the brightest stars. The reason seems to be 

not that Sandage's function differs somewhat from Salpeter' 

original function, but that the evolutionary theory used 

to obtain Y. , differs for Sandage's (and Salpeter ’s) 
bol 

function and for Idlis's rederived Salpeter function. 

Idlis also compares his Y^ ^ with observations of 
the Cygnus and Orion associations. These comparisons 
and others made by Usher with data from Trumpler's un- 
published catalogue of galactic clusters give reasonably 
good, agreement with the Idlis function. Walker's (1956, 
1957, 1961) extensive work on extremely young clusters 
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offers, in particular, further evidence of the secondary 
maximum. Van den Bergh's (1067) luminosity function does 
not show the secondary maximum, but he smoothed his data 
considerably and , investigated mainiy clusters which con- 
tain stars later than B5 and which show some evidence of 
evolution. The raw observed data of van den Bergh and 
of Sandage, however, do show a secondary maximum at 
“v • - 2 -° to - 2.5. Finally, the work of these two 
authors demonstrates the universality of the initial 
luminosity function among existing galactic clusters. 

7. Evolution of ini tial luminosity, f unction. — The 
consequences of evolving the initial Idlis function with 
the data given in Table 1 may be seen on Figure 1. (It 

was assumed that stars disappear from the main sequence 

18 

when X < 0.1.) Results are shown for values of k x 10 

18 

equal to 0, 0.22, and 1.0 gm/erg, and for k x 10 /L^ ^ 

2 

equal to 1.0 gm-sec/erg , where I* 13 is the initial 
luminosity of a star of 13 

The first case corresponds to homogeneous evolution 
without mass loss. The initial rapid brightening of the 
more massive stars is evident as hydrogen becomes quickly 
consumed. Their subsequent disappearance from the main 


12 


sequence occurs simultaneously, with .the. slow brightening 

• ' a 

of the less massive stars. As a result , after 3 x 10 

years a central accumulation of stars appears around 

Vi - - 3 - 


The choice of k = 0.22 x 10 


-18 


for the second case 


was directed by a desire to cause an initial decrease in 

\ 

L and then an increase back to L at half the initial 

o 

hydrogen content (cf. equation (26)). However, it should 

I : ■ ' 

be noted that this will not occur on the same time scale 
for stars of different masses. We have adopted Fesenkov's 
law (equation (1)) for the form of the mass loss rate. 

From the first case without mass loss, the most massive 

7 

stars begin to disappear at 2.5 x 10 years. Mass loss 

7 

prolongs their lifetime up to 5 x 10 years, but during 
this time they will have dimmed and brightened successively. 
Consequently the luminosity function shows no apparent 
change until their hydrogen is completely exhausted. 

Even smaller amounts of dimming and brightening account 
for the stationary behavior of the low-mass stars. Be- 
cause of the coarseness of the adopted it was im- 
possible to find a value of k that would make actually 

# • 

N show an initial decrease, subsequent increase, and final 


13 



decrease at high luminosities. 

-18 

Idlis's (1967) semi-mflpirieSkl value of k - 1 x 10 
corresponds to a mass loss which will dominate the evolu- 
tion. Adopting Fesenkov's law for the form of the loss 
rate, we obtain the results for the third case in Figure 1* 

On account of such an extensive loss, hydrogen depletion. 

\ ’ 

\ 

in the most massive stars is only fifty per cent even 
8 

after 5 x 10 years. The result is that the luminosity 

I . 

function loses weight continually at the brighter magni- 
tudes, while the secondary maximum gains and moves to 

7 

fainter magnitudes. By 5 x 10 years the accumulation 

toward = 0 has eliminated the secondary maximum 

8 

entirely. A new minimum appearing after 1 x 10 years 
is caused by the accelerating evolution of stars initially 
at * 0. These results are in many ways similar to 

those derived by Usher (1963) using only equations (1) 
and (2) . The adoption of varying X. and v, however, does 
cause the temporary disappearance of the central minimum 
in the luminosity function. 

The fourth case in Figure 1 illustrates the choice 
of a mass loss rate proportional to L^. The numerical 
coefficient corresponds to that used in the previous 
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case but is normalized to the initial luminosity of the 

most massive stars (13 Ih^) . The offset is to deplete 

significantly only the mass of the more massive stars 

8 

and hence to lower their luminosities. By 5 x 10 years 

* 

their hydrogen has been exhausted, but the onset of hydro- 
gen depletion in the stars of lower mass begins to in- 
crease these luminosities somewhat; a shift of the maximum 

= 0 to M bol a - 2 

, Idlis's rate of mass loss corresponds to about . 

—6 

10 lU^/year for the most massive stars evolving at their 

maximum (initial) luminosity. Underhill's (1960) value 
—7 

of 10 IftVyear derived from observations of Be shell 
stars is probably a very generous estimate. Hence Idlis's 
rate may be taken as an upper limit. 

8. .Comparison with o bserved luminosity functions . — 

Data for four galactic clusters (h Persei, the Pleiades, 

, Praesepe, and the Hyades) were taken from Sandage (1957a). 

For consistency with the adopted bolometric correc- 
tions were based on Idlis's work, as follows: from his 

f < 

Figure 1 the relation between spectrum and may be . 

read; then Parenago's relation between spectrum and B.C. 


' ’• . ‘ Vy ; 


in the luminosity function from 
is the result. 
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(Zonn & Rudnicki 1959) gives us the required My-B.C. re- 
lation. For My < - 4# the bolometric correction was as- 
sumed to be - 3. 

Figure 2 shoes normalized luminosity functions for 
the four galactic clusters down to My w + 2, with the 
omission of two extremely bright stars in h Persei. 
Comparison of the luminosity function of h Persei with 
Figure 1 shows reasonable agreement with almost any 
slightly evolved theoretical function. Indeed, even on 
the basis of inhomogeneous evolution Hayashi & Cameron 

7 

(1962) derive the young age of < 2 x 10 years. However, 
the tertiary maximum around = - 7 suggest inhomo- 

geneous evolution for at least some of the stars. 

The luminosity function of the Pleiades may be 

8 

reasonably well reproduced by the cases 0.5 - 1 x 10 

7 

years for k = 0 (no mass loss) or by the case 1 x 10 

—18 

years for k = 10 (maximum mass loss) . The former 
age is, of course, close to that given directly by the 
assumption of inhomogeneous evolution (Sandage 1957b) , 
since stars do not move far from the main sequence until 
the central hydrogen content is low, anyway. 

Praesepe and the Hyades show considerable evidence 


16 


“18 

of evolution# and on the mass loss theory (k ** 10 ) 

ima«« have •« ■ *. 10 7 and i m 10® years# respectively. 

However# in order for all the massive stars to have 
evolved away and for stars of lower mass to have brightened 

a 

to = 0, the ages of these two clusters must be in 

excess of 5 x 10 years in the case of no mass loss 
, (k • 0). Bandage (1957b) gives ~ 1 x 10 9 years for 
these clusters on the inhomogeneous theory. 

We conclude that cluster luminosity functions per ££. 
may not be used to decide for or against the theory of 
homogeneous evolution with mass loss# since any observed 
function may be reasonably well reproduced on this theory. 
Moreover# it does not seem possible to determine k and 
the age uniquely# by using such functions. On the other 
hand# the theory of inhomogeneous evolution without mass 
loss makes definite predictions of the luminosity function. 
Good evidence exists that this theory predicts correctly 
(Hayashi & Cameron 1962). Further observational arguments 
may also be raised in its favor (Stothers 1963b) . 

Apart from an effectively negligible observed rate 

f 

of mass loss (Stothers 1963a), two other direct observa- 
tional tests may be made of the Russian theory. A Ih— |j, 


relation should exist among the members of a star cluster# 
and improved cluster expansion ages should yield defini- 
tive ages, which, then narrow considerably the choice of 
comparative theoretical luminosity functions. 
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TABLE 1 

Adopted initial values 

v 

6 
6 
6 
7 

4.5 3.5 7 

3.5 , 3.5 7 

2.5 4 9 

2.0 4 9 

1.7 4 9 





*bol 

0.04 

0.06 

0.06 

0.05 

0.09 

0.18 

0.30 

0.13 

0.00 
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Fig. 2. 
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Evolution of Idlis's initial luminosity 
function for the following values of the 
mass loss parameters: i =1/(1) k = 0, 

(2Xfc - 0.22 x IQ -18 , (3) k * 1 x 10* 18 ; 

■ \ 

and t*2, (4) k - 1 x 10“ 18 L 13 . Ordinata 
intervals are equal to 0.05. Time, t, is 
given in years. 1 ; 

Normalized luminosity functions of four 
galactic clusters. 








